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Abstract 

Finite-element micromagnetic simulations are employed to study the chiral symmetry break- 
ing of magnetic vortices, caused by the surface roughness of thin-film magnetic structures. An 
asymmetry between vortices with different core polarizations has been experimentally observed for 
square-shaped platelets. E.g., the threshold fields for vortex core switching were found to differ 
for core up and down. This asymmetry was however not expected for these symmetrically-shaped 
structures, where both core polarizations should behave symmetrically. Three-dimensional finite 
element simulations are employed to show that a small surface roughness can break the symme- 
try between vortex cores pointing up and down. A relatively small sample roughness is found 
sufficient to reproduce the experimentally observed asymmetries. It arises from the lack of mirror- 
symmetry of the rough thin-film structures, which causes vortices with different handedness to 
exhibit asymmetric dynamics. 



PACS numbers: 75.75.+a, 76.50.+g, 07.85.Tt 



I. INTRODUCTION 



A magnetic vortex is a typical ground state of soft magnetic thin-film platelets with 
dimensions in the micron or submicron range and, e.g., a square or circular shape. The 
magnetostatic field tries to align the magnetization parallel to the sample boundaries, so it 
is forced into the sample plane and curls around the structure-center, forming the vortex. 
At the center of the vortex, however, the exchange interaction causes the magnetization to 
turn out of the sample plane, avoiding an antiparallel alignment of neighboring moments. 

mm- 

The in-plane curling magnetization can have two possible senses of circulation c: clock- 
wise or counterclockwise. These states are denoted by c = — 1 and c = +1, respectively. 
Analogously, the out-of-plane magnetization of the vortex core can have two possible 
polarizations p: either up or down, denoted by p = +1 and p = —1, respectively. The 
combination of c and p also gives the vortex a three-dimensional handedness. A vortex with 
c X p = +1 is defined to be right-handed, and one with c x p = — 1 to be left-handed [3]. 

Magnetic vortices are interesting configurations from a fundamental point of view, but 
are also possible candidates for future memory applications. The core polarization p, which 
can store one bit of information, can be switched by low-amplitude magnetic fields [1]. 
In-plane alternating [3] or rotating [5] magnetic fields can switch the core easily when they 
resonantly excite the so-called vortex gyrotropic mode. This is the mode that corresponds 
to a circular movement of the vortex core around the structure center [HI Ej. It typically has 
a resonance frequency in the sub-GHz range and a width in the order of tens of MHz [S] . 

In a square or disk-shaped sample, all combinations of c and p are symmetric, i.e., 
these states have the same energy and symmetric dynamics. In particular their gyrotropic 
eigenfrequency and switching thresholds are exactly equal [9J. This symmetry follows from 
the mirror-symmetry of the circular or square sample shape. E.g., the p = +1 and p = — 1 
states can be transformed into each other by mirroring the sample along the sample plane 
(defined to be the xy-plane, z is the out-of plane direction). Similarly the c = +1 and 
c = — 1 states can be transformed into each other by mirroring along the xz or yz plane. 
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However, several experiments on Permalloy (Ni8oFe2o) square platelets have revealed 
a relatively large asymmetry in the dynamics of a vortex core when its polarization was 
inverted. First of all, during a steady-state gyration of the vortex core, induced by an 
in-plane alternating field, it was found that the gyration radius was significantly different 
for the two core polarizations (TU]. Secondly, the required amplitudes for switching with 
rotating fields were found to differ for the two polarizations: the switching thresholds of the 
investigated sample differed by more than 10% for p = ±1 [5]. 

Since the investigated structures were geometrically symmetric, a possible explanation 
for the p-symmetry-breaking was sample roughness, which breaks the mirror- symmetry on 
a small length scale. This was suggested as a possible explanation in the abovementioned 
papers [HI HO], although it remained unclear if a small roughness could cause such a large 
p-asymmetry. 

Surface roughness is usually inevitable in evaporated or sputter-deposited thin films, so 
its influence on the vortex dynamics may have to be taken into account for the development 
of, e.g., vortex-based memories [Hj. Especially the different vortex core switching thresholds 
for core up/down [5] may have implications for such devices. 

II. MICROMAGNETIC SIMULATION SETUP 

In this work, a micromagnetic study of symmetry breaking in vortices by a chiral 
sample geometry is presented, with surface roughness as a particular case. An in-house 
developed 3D finite-element program [15] was employed to solve the Landau-Lifshitz 
equation [T2]. The finite elements allow the mesh to be irregular, so that arbitrary sample 
shapes can be well approximated. In order to efficiently calculate the magnetostatic field 
on such a mesh, the Fast Multipole Method [13] was used. Special care was taken so 
that the small truncation errors of this method would not break the mirror-symmetry of 
the sample. Material parameters typical for Permalloy were chosen: exchange constant 
A=13 X 10^^^ J/m, saturation magnetization Mj5=736 x 10^ A/m [H], damping parameter 
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a=0.01, anisotropy constant Ki=0, gyromagnetic ratio 7=2.211 x 10^ m/A s. A 3D 
discretization was performed with a cell size of approximately 3.9 nm x 3.9 nm in x and y, 
and four layers in z, corresponding on average to 12.5 nm (see Fig. |2]^c)). 

III. SIMULATION OF CHIRAL GEOMETRIES 

To understand the properties of magnetic vortices in chiral nanostructures, a first series 
of simulations has been carried out on three deformed platelet geometries, shown in row 1 
of Fig. [l] These were 250 nm x 250 nm wide and 50 nm thick. Sample (a) is a screw-like 
structure which is not mirror-symmetric, and we arbitrarily define its chirality to be 
left-handed. Sample (b) is obtained by mirroring sample (a) and is therefore defined as a 
right-handed structure. Sample (c), on the other hand, also has two "jumps", similar to (a) 
and (b), but this sample is mirror-symmetric. 

The gyrotropic eigenfrequencies and ground-state energies of these samples were 
simulated for the four possible combinations of c and p, and are listed in Tab. [Tj In the 
mirror-symmetric sample (c) all combinations of c and p yielded the same properties, 
while a c- and p-asymmetry was found for the chiral samples (a) and (b). In those cases, 
the eigenfrequency and ground-state energy changed when the vortex core polarization 
p was inverted. The exact same asymmetry was observed between states with different 
circulations c, and the asymmetry vanishesd again when both c and p were inverted at 
the same time. This residual cp-symmetry is fundamental and follows immediately from 
micromagnetic theory: inverting both c and p corresponds to inverting the magnetization 
m, which in the absence of an external field does not change the energetics of a magnet. 
The c and p asymmetries can thus be interpreted in terms of chirality: vortices with equal 
handedness c x p are fundamentally symmetric, but in a chiral sample geometry, vortices 
with opposite handedness are asymmetric. 

Note that mirroring the handedness of the vortex has the same effect on its properties 
as mirroring the geometry of the sample. E.g., a left-handed vortex in a left-handed 
sample (a) has the same properties as a right-handed vortex in a right-handed sample (b). 



4 



When we denote a left-handed structure by /i = +1, a right-handed structure hy h = —1 
and a mirror-symmetric structure hj h = 0, then all configurations with equal c x p x h 
are symmetric. In particular, mirror-symmetric structures are always symmetric for all 
combinations of c and p, and in all cases the fundamental cp-symmetry exists: vortices with 
equal handedness are always symmetric. 

In the absence of a magnetic field, a chiral sample geometry was thus found to be 
necessary for breaking the c and p-symmetries. This necessary condition can be intuitively 
understood by investigating out-of-plane component of the magnetization. This is shown 
in Fig. [T| for the three respective sample geometries. Due to the shape anisotropy, which 
wants to align the magnetization parallel to the sample boundaries, the magnetization 
acquires an out-of-plane component 7^ in the regions where the films are deformed. 
These regions couple to the vortex core via the exchange interaction and the dipole-dipole 
coupling. As a consequence, vortex cores with different polarizations can have different 
energies in such an environment. When either the in-plane circulation c, or the core 
polarization p is inverted, the coupling changes sign as well. The asymmetry can thus be 
understood as a cp-coupling due the chiral sample geometry. 

Row 2 of Fig. [1] shows that in sample (a), the left-handed circulation causes the extra 
out-of-plane magnetization at each of the jumps to point up. The coupling between those 
regions and the vortex core happens to be favorable in this case, causing this configuration 
to have the lowest energy. In the right-handed structure (b), the opposite takes place. 
The out of-plane magnetization at the jumps is antiparallel to the vortex core, causing the 
configuration to have a higher energy. On the other hand, if the sample is mirror-symmetric 
like sample (c). The magnetization points up at one jump and down at the other one, so 
that the coupling is equal for core up and down. 

Note that the presence of an external field in the sample plane can cause an additional 
symmetry breaking. E.g., in the mirror-symmetric sample (c), a static field in the 
x-direction will displace the vortex core down, towards the region with a positive out-of 
plane magnetization. It can easily be understood that an up vortex will behave differently 
from a down vortex in such an environment. In this particular case, the energy of the up 
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vortex was found to be lower. However, even in the presence of an in-plane magnetic field, 
the deformed geometry is still necessary to break the symmetry between the four vortex 
states. 

IV. SIMULATION OF ROUGHNESS 

The previous set of simulations was artificially chosen to illustrate the mechanism behind 
the c and p symmetry-breaking by a lack of geometric mirror-symmetry. Random surface 
roughness is a particular case where the mirror-symmetry is (almost always) broken, and 
where the condition for a cj9-asymmetry is thus fullfilled. How large this asymmetry is and 
which vortex handedness is favored, can be expected to depend the details of the roughness. 

To illustrate that a chiral sample geometry can explain the experimentally observed 
asymmetries, we performed a second set of calculations of samples with a random surface 
roughness, similar to the samples investigated in the abovementioned experiments [TU] . 
The simulated geometry was inspired on a scanning electron micrograph of a typical 
Permalloy sample investigated in those experiments, shown in Fig. |2]^a). The sample lies on 
top of a relatively rough copper stripline, which is used for the magnetic excitation [H |5l [10] . 
Therefore, a square platelet was simulated as if it were deposited on a rough underground 
while conserving its thickness. A peak amplitude of 10 nm was chosen for this roughness. 
To make the geometry more realistic, a smaller additional roughness of 5 nm was added to 
the top surface, although this did not significantly change the outcome of the simulations. 
The length scale of the roughness was also chosen to mimic the electron micrograph in Fig. |2} 

Due to the random nature of the roughness, the effect in different samples can be 
expected to have a different strength or might incidentally even vanish. It lies however 
not in the scope of this work to study the differences between many rough geometries, 
neither to exactly reproduce the detailed roughness of the experimentally investigated 
structures. We only want to illustrate that realistic roughness parameters are sufficient 
to explain the experimentally observed asymmetries. Therefore, we have simulated three 
rough geometries based on different random number seeds, and we present the simulation 
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with the median symmetry breaking between the gyrotropic eigenfrequencies in Fig. |3} The 
expected asymmetry between the states with different vortex handednes c x p was found. 
The figure also shows the out-of plane moments that arise due to the roughness. A side 
view of a vortex structure in a rough thin film (Fig. |3] (a)) illustrates how these out-of-plane 
moments are especially pronounced when the full layer happens to be bent upwards or 
downwards, which is exactly the case when the film is deposited on a rough underground. 
The expected asymmetry was also found in the ground-state energy: for c x p = +1 this 
was 3.0151x10-^^ J, while for c x p = -1 it was 3.0102x 10"^^ J. 

The observed difference in resonance frequencies of about 20 MHz is in the order of the 
width of the resonance, and is thus expected to significantly affect the dynamics. E.g., 
when the excitation is chosen in resonance for the up vortex, it may already be well out 
of resonance for the down vortex. The small asymmetry in gyration frequency could thus 
cause a relatively large asymmetry in the response to an alternating field near the resonance 
frequency. 

This has been quantified by applying a rotating magnetic field in gyrotropic sense of 
the vortex core [5J to each of the c,p states. The field was slowly ramped up with a time 
constant of 15 ns, until it was strong enough to switch the vortex core polarization. A 
frequency of 908 MHz was chosen, which lies closest to the eigenfrequency of the cp = +1 
configurations. For the different combinations of c and p, switching thresholds between 
4.19mT and 4.69 mT were found. This asymmetry between the switching thresholds is of 
the same order of magnitude as the experimentally observed asymmetry of about 10% [3]. 
The presence of a small roughness can thus explain the experimental asymmetries. Note 
that in the absence of roughness, all combinations of c and p yielded exactly the same 
switching threshold and gyrotropic eigenfrequency. 

Finally, it was asserted that the sample size has no critical effect on the symmetry 
breaking, by simulating three 500nm x SOOnmx 40nm squares with random roughness. 
Due to computational limitations, these could only be discretized in two dimensions, but 
an asymmetry of the same order of magnitude was found. This reflects the local nature of 
the effect of the roughnss on the vortex core. 
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V. CONCLUSION 



We have presented micromagnetic simulations showing that the absence of mirror- 
symmetry in the sample shape can break the symmetry between vortices with opposite 
handedness c x p. More specifically, inverting either the in-plane circulation c, the core 
polarization p or the chirality of the sample shape changes the vortex energy and dynamics 
in an asymmetric way. A particular case of chiral sample geometry is a random surface 
roughness. The out-of-plane moments caused by the roughness couple with the vortex core, 
breaking the symmetry between core up and down. Although the global symmetry breaking 
by the roughness is small, the coupling of the vortex core with the local out-of-plane mo- 
ments is strong enough to explain the experimentally observed asymmetries. Although we 
have shown that roughness can break the chiral symmetry of vortices, this work does not 
exclude the possibility of other, additional sources of symmetry breaking. 
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(a) left-handed (b) right-handed (c) mirror-symmetric 




FIG. 1: Row 1: deformed thin-film platelets. The structures (a) and (b) are chiral and are defined 
to be left-handed and right-handed, respectively, while sample (c) is mirror-symmetric and thus 
has no handedness. Row 2: out-of-plane component of the magnetization in those structures, 
for c = +1, = +1. In the chiral structures (a) and (b), niz at the position of the jumps points 
either twice up or twice down, so no symmetry exists between the magnetization of (a) and (b). 
In the mirror-symmetric structure (c), niz points once up and once down, so that here a symmetry 
between the two possible p-states remains. 




FIG. 2: (a) Scanning electron micrograph of a typical 500 nm Permalloy platelet in which a p- 
asymmetry was observed. The platelet lies on top of a rough copper stripline. (b) Simulation 
geometry of a 250 nm platelet with a roughness which mimics the effect of a rough underground, 
(c) Detail of (b), showing the 3D finite-element discretization. 



10 




FIG. 3: (a) Out-of-plane component ruz of the magnetization in a rough sample, for the possible 
combinations of c and p. Additional out-of-plane moments arise, which are not present in perfectly 
flat structures. Their direction is dependend on the circulation c. These out-of-plane moments 
interact with the vortex core via the exchange and dipole-dipole coupling, resulting in different 
energies for a vortex that points up or down, (b) Illustration how these out-of-plane moments 
arise. The shape anisotropy keeps the magnetization parallel to the sample edges, locally tilting 
the magnetization up or down at some positions. The resulting out-of-plane moments are especially 
large when the top and bottom surface are nearly parallel, as is caused by a rough underlayer. 

Tables 
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sample (a): left-handed geometry 



cp 


1 


-1 


fo 


914.9 MHz 


837.8 MHz 


Eo 


2.97x10-1^ J 


2.99x10"^^ J 




cp 


1 


-1 


fo 


837.8 MHz 


914.9 MHz 


Eo 


2.99x10-1^ J 


2.97x10-1^ J 




cp 


1 


-1 


fo 


1054.4 MHz 


1054.4 MHz 


Eo 


2.94x10-1^ J 


2.94x10-1^ J 



.e (b): right-handed geometry 



sample (c): mirror-symmetric 



TABLE I: Simulated resonance frequencies fo of the vortex gyrotropic mode and ground state 
energies Eq, for the 250 nm x 250 nm x 50 nm structures shown in Fig. [T] The four combinations 
of vortex circulation c and polarization p were simulated; states with equal cx p yielded the same 
results. Structures (a) and (b) are chiral and are each other's mirror image, while structure (c) is 
mirror-symmetric (see Fig. [T]) . An asymmetry between vortices with a different handedness c x p 
arises only for the chiral structures. 
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